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Abstract 

We compute the electromagnetic form factor of the pion using non-perturbatively 
O(a) improved Wilson fermions. The calculations are done for a wide range of pion 
masses and lattice spacings. We check for finite size effects by repeating some of the 
measurements on smaller lattices. The large number of lattice parameters we use allows 
us to extrapolate to the physical point. For the square of the charge radius we find 
(r 2 ) = 0.441(19) fm 2 , in good agreement with experiment. 



1 Introduction 

For some time now it has been possible to explore the structure of hadrons from first principles 
using lattice QCD. Since the pion is the lightest QCD bound state and plays a central role in 
chiral symmetry breaking and in low-energy dynamics, a thorough investigation of its internal 
structure in terms of quark and gluon degrees of freedom should be particularly interesting. 
We have started to explore the structure of the pion in a framework using generalised parton 
distributions, or more precisely their moments pQ. As a generalisation of parton distributions 
and form factors they contain both as limiting cases. In this work we restrict ourselves to 
results for the pion electromagnetic form factor from Nf = 2 lattice QCD simulations, 
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based on 0(a) improved Wilson fermions and Wilson glue. Initial studies on the pion form 
factor by Martinelli et al. and Draper et al. [2[ [3] were followed by recent simulations in 
quenched [H [7] and unquenched QCD [5] ■ In this work, we improve upon previous 
calculations by extracting the pion form factor for a much larger number of (3, k combinations, 
which allows us to study both the chiral and the continuum limit. Furthermore, two finite 
size runs make estimates of the volume effect possible. 

2 The pion form factor in lattice QCD 

The pion electromagnetic form factor F n describes how the vector current 

V^x) = | u(x)7 M u(x) - \2(x)^,j,d{x) (1) 

couples to the pion. Writing p and p' for the incoming and outgoing momenta of the pion, it 
is defined by 

(n+(p')\ V^O) \n + (p)) = + F 1T (Q 2 ) , (2) 

where the momentum transfer is = (jp' — Pn) and its invariant square is q 2 = —Q 2 . 

For our lattice calculation we want to simplify the flavour structure of Eq. ([T]) . Invoking 
isospin symmetry one finds 

(•7r + | |u7 m m - ^dl^d = (7T+| u-f^u \tt + ) = - dj li d\ir + ) . (3) 

It is hence sufficient to limit the calculation to a single quark flavour in the vector operator. 
We use the unimproved local vector current on the lattice; the corrections due to the improve- 
ment term [10] are quite small and will be discussed later. Since this current is not conserved, 
renormalisation has to be taken into account. Because in the forward limit (Q 2 = 0) the 
form factor is simply the electric charge of the pion, we can normalise our data appropriately. 
We can also use the known renormalisation constant Zy (taken for example from [llj ) as a 
cross-check for our simulation. 

To compute the matrix elements in Eq. ^ on the lattice, one has to evaluate pion three- 
point and two-point functions. We then apply a standard procedure to extract the pion form 
factor , where one constructs an appropriate ratio for the observable [12j [13] . Let us start 
by looking at the three-point function. The general form is given by the correlation function 

Capt(t,p',p) = (^(tsinkiP') u{thnu(t) r)l(t SOUIce ,p)) (4) 

and depicted in Fig. [TJ Here we denote the sink and source operators for a pion with given 
momentum and at given time-slice by 7? T (i s i n k,p') and r]l(t sourcB ,p), respectively. Using the 
transfer matrix formalism and inserting complete sets of energy eigenstates, the three-point 
function is then of the form 

r> n s*/ x\ I wnirfnl cnM tzw (°\^(p ') K(p')) (^(p^vUp) 1°) 
C 3 pt(t,p ,p) = {Tr(p )|u(0)7 m u(0) \tt(p)} 



x 



f e -E p ,(tsink-t)-E p t + (_ 1 )"4 c -E p ,(t-t sink )-E p (T-t)^ _j ^ ^ 

where T is the time extent of our lattice and 

n = I 1 f ° r ^ = 4 ' (61 

ni | otherwise. ^ ' 
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Figure 1: A sketch of the three-point function with the pion source at time 0, pion sink at 
^sink, and the operator acting at time t. 



Note that we have omitted excited states in Eq. ([5]) and already inserted our choice for the 
time-slice of the pion source, t aomcc = 0. We choose the sink of the three-point function as 
^sink = T/2, so that the correlation function is symmetric or antisymmetric with respect to 
this time, 

C 3pt (t,p',p) = (-1)"* C 3pt (T - t,p',p). (7) 

We can then separate the correlation function into contributions from t to the left and to the 
right of igink (referred to as l.h.s. and r.h.s. in the following) and neglect either the second 
or first term in Eq. ([5|) since it is exponentially suppressed in the regions of t from which we 
will extract the form factor. 

The two-point function has the form 

c 2pt (t,f) - (Q| rf) im^m m \% -^ 2coMEp{T/2 _*)]+..., (8) 

where again we omitted higher energy states. Comparing the two- and three-point func- 
tions ([5]) and ((5]), a ratio can be constructed that eliminates the overlap factors such as 
(0 ijwip') \ n (p')) an d partially cancels the exponential time behaviour appearing in Eq. 
This technique also has the advantage that fluctuations of the correlation functions tend to 
cancel in the ratio and we thus obtain a better signal. With our choice i s ink = T/2, such a 
ratio is 



R(t) C ^P'>P) 



C 2 pt(isink,P') 



C 2 pt(Aink ~ t,p) C 2pt {t,p') C2ptfaink,p / ) 



C2pt(isink - t,p') C 2pt (t,p) C 2 pt{t s ink,P) 



(9) 



Similar ratios have already been used in earlier works on pion and nucleon structure. Here we 
take the somewhat more complicated ratio (O , which was used for the nucleon in [12] , because 
we use momentum combinations with \p\ ^ \p'\- Contributions to this ratio from excited 
states with energy E' are suppressed as long as i s ink — t ^> l/(£" — E) and t 3> l/(£" — E) 
where E is the pion energy. A potential problem is that, due to the exponential decay of 
the pion two-point function, the signal at t = < s i n k for non-vanishing momenta is poor. For 
finite statistics the two-point function can then take negative values, which prevents one 
from evaluating the square root. We try to overcome this difficulty by shifting the two-point 
functions C 2p t(t,p) that enter with t = i s ink- Using the identity 

/, -k C*2pt (isink — ^shift , P) , in N 
2 pt(tsink,Pj = w„ > (lUj 

cosh(^pUift) 

valid for t S i n k = T/2, we shift by £ s hift = 6, which significantly reduces the number of 
negative two-point functions. Nevertheless there are still momentum transfers Q 2 for which 
the argument of the square root in the ratio © is negative. Those values are discarded when 
we evaluate the form factor. 
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Figure 2: Examples of the ratio R(t) in Eq. © on the 24 3 x 48 lattice at {3 = 5.25, 
K = 0.13575, multiplied with an appropriate sign factor (— l)" 4 for t > i S mk- The left plot 
shows a proper plateau in the forward case Q 2 = and the right plot shows the ratio for 
Q 2 = 0.31 GeV 2 where no plateau is expected (see text). The dashed lines indicate the 
regions we average over. 



For Q 2 ^ the ratio ((9|) docs not exhibit a proper plateau that could immediately be 
used for fitting. This is due to our choice for t S i n k, for which the time dependence of the pion 
two-point function cannot be approximated by a single exponential in the t regions we use 
to extract the form factor, see Eq. (jSJ). In fact, we now show that the ratio is approximately 
antisymmetric around the central point t = i S i n k/2 = T/4 of the l.h.s. (as well as around 
t = 3T/4 on the r.h.s.). Defining 5 = t — i S i n k/2 and expanding the ratio and its exponentials 
in Eq. (J9} around 6 = we find 

R(t) = C(E p ,E p/ ,Q 2 )[l + 25c s (E p ,E p ,) + 28 2 c 2 {E p ,E p/ ) + 0(8 3 )] , (11) 

where 



cs(E p ,E pl ) 



E p i E p 



X -)- Q E P Wk \ _|_ Q E p> *sink 



C(E P ,E P ,,Q 2 ) = { A±PA Fw{Q 2y (12) 

4y<E p ,E p 

When averaging R(t) in a symmetric interval around t = T/4, the antisymmetric piece 
proportional to c$ in pip drops out. However, such an averaged signal also includes unwanted 
symmetric contributions. Fortunately, for our pion masses and lattice momenta already the 
leading symmetric term is negligible, because with the lattice spacing a we have c 2 ~ 10~ 4 a~ 2 
and S 2 < 4a 2 in our fits. We hence obtain a good signal for the averaged ratio. The same is 
true for the r.h.s. ratio and its central point t = 3T/4. A typical ratio at non-zero momentum 
transfer is shown in Fig. [5] for one of our data sets, along with the familiar plateau for zero 
momentum transfer. Note that the ratio ^ does not exhibit a plateau for arbitrary momenta. 
To visualise the absence of possible contributions from excited states one has to consider the 
ratio for \p\ = \p'\. In this case the time dependence of the three-point function should 
vanish. We have checked that this is indeed the case in the region we average over, within 
the expected increase of noise for higher momenta or lower pion masses. 

From Eqs. (fTTj) and (fT2"|) we see that the lattice ratio @ can be used to extract the 
form factor F 7T (Q 2 ). Using then several combinations of momenta p and p' that all give the 
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Table 1: Overview of our lattice parameters. For physical units the Sommcr parameter with 
rp = 0.467 fm has been used. The error on m w is purely statistical. 
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4A 
jr 


K 


N 3 xT 


77} fGeVl 


d [fin] 


L [fm] 


1 y tra] 


5.20 


1 


0.13420 


16 3 x 32 


1.007(2) 


0.115 


1.8 


O(5000) 




2 


0.13500 


16 3 x 32 


0.833(3) 


0.098 


1.6 


0(8000) 




3 


0.13550 


16 3 x 32 


0.619(3) 


0.093 


1.5 


C(8000) 


5.25 


4 


0.13460 


16 3 x 32 


0.987(2) 


0.099 


1.6 


0(5800) 




5 


0.13520 


16 3 x 32 


0.829(3) 


0.091 


1.5 


O(8000) 




b 


0.135 lb 


1 A 3 w AO 


0.59^(1) 


n no A 

0.084 


2.0 


/ C A AA\ 

C(5900) 


5.26 


7 


0.13450 


16 3 x 32 


1.011(3) 


0.099 


1.6 


0(4000) 


5.29 


8 


0.13400 


16 3 x 32 


1.173(2) 


0.097 


1.6 


0(4000) 




9 


0.13500 


16 3 x 32 


0.929(2) 


0.089 


1.4 


0(5600) 




10 


0.13550 


24 3 x 48 


0.769(2) 


0.084 


2.0 


O(2000) 




11 


0.13590 


24 3 x 48 


0.591(2) 


0.080 


1.9 


0(5900) 




12 


0.13620 


24 3 x 48 


0.400(3) 


0.077 


1.9 


0(5600) 


5.40 


13 


0.13500 


24 3 x 48 


1.037(1) 


0.077 


1.8 


0(3700) 




14 


0.13560 


24 3 x 48 


0.842(2) 


0.073 


1.8 


0(3500) 




15 


0.13610 


24 3 x 48 


0.626(2) 


0.070 


1.7 


0(3900) 



same Q 2 provides an over-constrained set of equations, from which we determine F V {Q 2 ) 
by \ 2 minimisation. We increase the quality of our signal by averaging the ratio over the 
contributions on the l.h.s. and r.h.s. This requires the additional sign factor (— l)™ 4 between 
the two sides, as can be seen in Eq. ([7J)- The energies E p and E p > appearing in (fTQ|) and (jXSJl 
are calculated using the lattice pion masses and the continuum dispersion relation. We also 
performed a test of the dispersion relation for some of our lattices. It was increasingly difficult 
to extract a signal for higher momenta, especially for the lowest pion masses. However, we 
found that the continuum dispersion relation can be used to describe the data and that a 
lattice dispersion relation is not favoured. 

3 Simulation details 

We perform our simulations with two flavours of non-perturbatively clover-improved dy- 
namical Wilson fermions and Wilson glue. Using these actions, the QCDSF and UKQCD 
collaborations have generated gauge field configurations with the parameters given in Ta- 
bic [1] where we have used the Sommcr parameter with tq — 0.467 fm (see |14j and |15j ) to 
set the physical scale. This large set of lattices enables us to extrapolate to the chiral and 
the continuum limit. For two sets of parameters (/? = 5.29, k = 0.1355,0.1359) we also have 
a choice of lattice volumes (12 3 x 32, 16 3 x 32 and 24 3 x 48) in order to study finite volume 
effects. 

Starting with the lattice version of the three-point function, Eq. (U), we follow [TB] and 
find that it is sufficient to calculate 

£ ^e-^V^Tr TG(y,z) 1 ,G(z 1 x)^G{x,y)) g (13) 
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with X4 = 0, 2/4 = T/2, z± = t. Here G(y, z) is the fermion propagator, the average is taken 
over the gauge fields, and the trace is over the suppressed Dirac and colour indices. The 
matrix T represents the Dirac structure of the pion interpolating field T)„, while the Fourier 
transformations ensure that we have fixed momenta at the operator insertion and the sink. 

We use two different pion interpolating fields to create the pions on the lattice, namely 
a pseudo-scalar and the fourth component of the axial- vector current, which both have the 
correct quantum numbers. For a given momentum p they read 

rfe (t, p ) = c~' lp ' S d(x) T u(x) , r = 75 or 7475 (14) 

X 

with X4 = t. We apply Jacobi smearing [17] at the source as well as the sink to increase the 
overlap of the lattice interpolating fields with the physical pion states. 

The three-point function (fT5|) is then evaluated by applying the sequential source tech- 
nique as indicated in Fig. [TJ This makes it efficient to use a large number of momentum 
transfers, as required for calculating form factors. A large set of momenta is necessary to 
assess the Q 2 dependence, and having several combinations of p' and q belonging to the same 
Q 2 makes the fits more reliable. We use three final momenta p' and 17 momentum transfers 
q, giving a total of 51 combinations for an over-constrained fit for F n at 17 different values 
of Q 2 . In units of 2ir / ' L the momenta are given by 

f = (0,0,0), (0,1,0), (1,0,0), 

q = (0,0,0), (-1,0,0), (-1,-1,0), (-1,-1,-1), (15) 
(-2, 0,0), (-2,-1,-1), (-2, -2,-1), •■• 

where • • • stands for all possible permutations w.r.t. the components. The errors we quote 
for our results are statistical errors obtained by the jackknife method. 

4 Experimental data for the pion form factor 

Let us now take a brief look at the experimental measurements of F. K {Q 2 ) to which we compare 
our lattice results. Very accurate data up to Q 2 = 0.253 GeV 2 have been obtained in [18] 
from elastic scattering of a pion beam on the shell electrons of the target material. At higher 
Q 2 the pion form factor has been extracted from ep — ► enir + , which is considerably more 
involved (see [THj for a recent discussion). We only use here data from [SO] [2T] [22] , where the 
cross sections for longitudinal and transverse photons have been experimentally separated 
by the Rosenbluth method^ Together these data span a range from Q 2 = 0.35 GeV 2 to 
2.45 GeV 2 . 

We find the experimental data on F v well described by a monopole form 

^(Q 2 ) = T^/^ (16) 

with a fit of the combined data from [H [501 123 H2 giving M = 0.714(4) GeV at x 2 /d.o.f. = 
1.27. This is remarkably close to the result M = 0.719(5) GeV at x 2 /d-0.f. = 1.13 obtained 
when fitting only the data of [TH] with its much smaller range in Q 2 , which illustrates the 
stability of a monopole form up to 2.45 GeV 2 . 

^^For the data from [21] we use the results of the re-analysis in |23| 
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Table 2: Values of the squared pion charge radius obtained from different data sets for 
F- K {Q 2 ) and with different fitting assumptions. Details of the fits are given in the text. 



data (r 2 ) [fm 2 ] 

global average, PDG 2004 [Mj 0.452(11) 

Amendolia [IS], fit 1 0.439(8) 

fit 2 0.431(10) 

fit 3 0.451(6) 

combined data [TJ H51 [HI H2] 0.458(5) 



The low-Q 2 behaviour of is characterised by the squared charge radius 



(r 2 ) 



d^(Q 2 



For a monopolc form ()16|) one has 



dQ 2 



(r 2 ) = 6/M 2 . 



Q 2 =0 



(17) 



(18) 



In Tabled] we list the values obtained from a number of fits to F^. The PDG average [H] 
uses results from form factor data at both spacelike and timelike virtualities. The three fits 
to the Amendolia data jTHJ illustrate that different fitting procedures can give results with a 
variation much bigger than the quoted statistical and systematic errors. Fit 1 (whose result 
is the one retained in the PDG average) is based on a representation of F„ as a dispersion 
integral. Fit 2 was also given in [18| and assumed a monopole form (|16[) with a normalisation 
factor allowed to deviate from 1 by ±0.9%, which corresponds to the overall normalisation 
uncertainty of the measurement. Fit 3 assumes a monopole form with normalisation fixed to 
1, as does the fit to the combined data of [TH ] 120 ] |2"T ] W\ . 



5 Results 

5.1 Fits to lattice data and extrapolation in m n 

We start the discussion of our results by explaining our fitting procedure, including combined 
fits to all data sets. In the next subsection we will argue that lattice artifacts are small. To 
obtain the physical form factor we have to renormalise our lattice result, F™ n ~ ZyF^ xc . 
As mentioned in Section [2] we can do this by using the electric charge of the pion as input, 
i.e. 

F^ at ' baro (0) 

so that i^. at < ron (0) = FP hys (0) = 1. We then use a monopole ansatz to fit the actual data for 
the form factoi] 2 ] 

F " ( « 2 > = TToW (20> 



z^lat,bare \ 

i7ilat,ron ( ' / ^ ) 2\ _ fjr vV I { -i n \ 



2 We will from now on use the renormaliscd values and drop the superscripts unless required. The supcr- 
and subscripts 'lat' and 'phys' respectively refer to observables at lattice pion masses and at the physical 
point. 



7 



where we have M\ ai as a fit parameter for each of our lattices at its lattice pion mass m^iat. 
The quality of this fitting ansatz will be discussed below. 

Using this fitting function, we compare the results obtained with the two pion interpo- 
lating fields (fT4"]) and observe several differences. In general, the matrix elements for pions 
using r = 7475 display a slightly cleaner signal with more data points in Q 2 , i.e. less con- 
tamination due to negative two-point functions. Fitting the monopole form (|20j) to the form 
factor for both pion interpolators we find that the x 2 /d.o.i. differs on average by about a 
factor of 2, ranging from 0.18 - 1.72 (0.23 - 3.49) for the interpolator with 7475 (75). The 
fitted monopole masses for the r = 75 pions lie consistently above the ones for T = 7475 but 
agree within errors for most lattices. In an exploratory extraction of the pion energies from 
the two-point functions with non-vanishing momentum on a sub-set of our lattices, we also 
found that the pseudo-scalars with T = 75 had a worse signal at higher momenta. A similar 
observation was made in [5] and may explain the difference in quality of the form factors 
extracted from the two pion currents. Because of the better signal, we will mainly discuss 
results for the pions created with T = 7475 in the remainder of this work. 

To obtain the pion form factor at the physical pion mass we extrapolate the values for 
M[ a t, given in Table [51 to the physical point. We tried different extrapolations in the square 
of the pion mass, see Table including also a fit inspired by chiral perturbation theory and 
used in [S] . For the latter we chose the fit range of m 2 lat < 0.8 GcV 2 . Varying this fit range 
within reasonable bounds did not have a significant effect on the extrapolated value of M p i iys . 
We find the best \ 2 value for fit 2, where M 2 at depends linearly on m 2 . The extrapolations in 
the remainder of this paper are based on this ansatz. We will however include an estimated 
systematic error of AA/ 0Xt = 35 MeV from the difference of fits 1 and 2 in our final result 
(this is bigger than the difference between fits 1 and 4, whereas fit 3 gives a significantly 
worse description of the data) . Figure [3J shows the extrapolation to the physical pion mass 
based on fits 2 and 4. We remark that our lattice with the lowest pion mass, m n = 400 MeV, 
is completely consistent and increases our confidence in the fit and fit ansatz. However, due 
to the larger statistical errors it has little weight in this result: when leaving it out of the 
fit Mphys changes only by 1 MeV. The corresponding run and several others at small pion 
masses arc still in progress. It is obvious that one needs higher statistics for this point to be 
significant. 

We include the m w dependence of the monopole mass of fit 2 in a combined fit to all our 
lattice data available. This fit has the same monopole form as in (|2"D)) with one additional 
parameter to incorporate the m w behaviour, 



The two fit parameters, cq and c\, describe the relation between the monopole mass and 
the pion mass, and we immediately obtain the form factor F^ hys (Q 2 ) = F v (Q 2 , m 2 , pw) m 
the physical limit. The fitted parameters are c$ = 0.517(23) GeV 2 and ci = 0.647(30) with 
X 2 /d.o.L = 0.64. This gives M phys = M(m 2 phys ) = 0.727(16) GeV, in good agreement with 
the experimental result. 

Figure U] shows experimental data along with the combined fit with its extrapolated 
curve. For this plot, our data at the lattice pion masses is shifted to the physical pion mass 
and plotted on-top of the extrapolation. We do this by subtracting from the individual 
lattice points, F\ &t (Q 2 ), a value (F n (Q 2 ,m 2 ut ) - F n (Q 2 ,m 2 phys )) calculated with the fit 
parameters of Eq. (f2"Tj) at the respective pion masses. The errors are left unchanged. We find 
good agreement between our simulation and the experimental results. This is emphasised by 




W ,mj- 1 + Q2/ M 2 (m 2) 
M 2 (ml)=c +c l ml. 



(21) 
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Table 3: Monopole masses M\ at obtained from fits to ([20]) for each of our lattices. The last 
column gives an estimate for the shift AMi at = M(m 2 , oo) — M(m%, L) of the monopole mass 
due to finite volume effects. It is obtained from the empirical fit ([25]) discussed in Section 15721 



It 


771 [GeVl 


L [fm] 




M, at [GeVl 


AMt at \MeV] 


1 


1.007(2) 


1.8 


9.4 


1.104(22) 


0.3 


2 


0.833(3) 


1.6 


6.6 


0.997(21) 


4.3 


3 


0.619(3) 


1.5 


4.7 


0.880(24) 


35.2 


4 


0.987(2) 


1.6 


7.9 


1.089(20) 


1.1 


5 


0.829(3) 


1.5 


6.1 


0.975(17) 


7.2 


6 


0.597(1) 


2.0 


6.1 


0.870(22) 


8.0 


7 


1.011(3) 


1.6 


8.1 


1.066(25) 


0.9 


8 


1.173(2) 


1.6 


9.2 


1.157(20) 


0.3 


9 


0.929(2) 


1.4 


6.7 


1.051(15) 


3.7 


10 


0.769(2) 


2.0 


7.8 


0.971(14) 


1.3 


11 


0.591(2) 


1.9 


5.7 


0.854(15) 


12.6 


12 


0.400(3) 


1.9 


3.8 


0.783(36) 




13 


1.037(1) 


1.8 


9.7 


1.099(13) 


0.2 


14 


0.842(2) 


1.8 


7.5 


0.981(14) 


1.8 


15 


0.626(2) 


1.7 


5.3 


0.847(17) 


18.2 



Table 4: Different forms used to extrapolate the monopole mass to the physical value of m^. 
In fit 4 we have L = l/(47r/) 2 log(m^ ut /fJ. 2 ), where \i = 1 GeV and f n « 92MeV is the pion 
decay constant. 



# 


extrapolation ansatz 


X 2 /d.o.f. 


Cl 


M phys [GeV] 


1 


Mat = c o + Ciml M 


1.31 


0.322(15) GoV" 1 


0.761(13) 


2 


M L = c o + cirn>l M 


0.93 


0.647(30) 


0.726(16) 


3 


l / M L = co + c x ml M 


3.25 


-0.575(31) GcV~ 4 


0.833(9) 


4 


®/ M L = co + c x ml Mt - L 


1.11 


-4.33(62) GcV~ 4 


0.715(4) 
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Figure 3: Extrapolations of the squared monopolc mass against the squared pion mass. The 
solid line with error band is linear extrapolation as obtained from fit 2 while the dotted line 
shows the central curve for fit 4 (whose fit range is limited to m 2 < 0.8 GeV 2 ). The cross 
marks the monopolc mass corresponding to the PDG value [24] of the pion charge radius, 
sec Eq. (fT8| and Table [2l The different symbols refer to our /3-values: squares (5.20), circles 
(5.25), half-full circle (5.26), diamonds (5.29), and hexagons (5.40). 



the insert in Fig. HI which shows the region Q 2 < 1 GeV 2 , where most of the experimental 
points lie. The same fit for the pions with T = 75 gives Af p h ys = 0.773(17) GeV, with a bigger 
X 2 /d.o.f. of 1.01. 

We now investigate the validity of the monopole ansatz for our data. Instead of constrain- 
ing the fitting function to a monopolc form, one can also take a general power law, i.e. use a 
function 

F 7r (Q 2 ,m 2 Hfl+ f ) 

\ pM 2 (ml)J (22) 

M 2 (m 2 ) = c + cito 2 , 

with an additional parameter, p. Note that the relation (TTS)) is still valid, independent 
of p. A combined fit to all our data sets results in p = 1.173(58), now with a mass 
-Mphys = 0.757(18) GeV and a x 2 /d.o.f. = 0.58, indicating that the monopole form is a 
good description. Taking the difference between this number and the result of the fit to (f2"Tj) . 
we assign a systematic error of AMj t = 30 MeV to M p h ys due to the ansatz for the fitting 
function. Another alternative is to calculate an effective monopole mass for every momentum 
Q 2 separately by solving Eq. (|20|) for M\ at : 

-1/2 

(23) 



Mcff(Q 2 ) = Q 



F^(Q 2 ) 
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0.5 1 1.5 2 2.5 3 3.5 4 



Q 2 [GeV 2 ] 

Figure 4: Combined fit to ([2"Tj) of our data for all lattices. We plot experimental data 
(diamonds) [THIHHHS] and lattice results extrapolated to the physical pion mass as explained 
in the text. To avoid having a cluttered plot we do not show lattice results with errors bigger 
than 80%, which are nevertheless included in the fit. The insert shows the good agreement 
to the experimental data for a momentum transfer of up to 1 GeV 2 . Also included is an error 
band for the fit. 

We show such effective masses for some of our lattices in Fig. [51 where one can see that the 
effective monopole masses stay constant within errors over a large range of Q 2 and agree with 
the monopole masses given in Table [3] This again indicates that the monopole is a good 
description for our data. The validity of the fit over the whole Q 2 range is further tested by 
combined fits to Eq. ([21]) in a limited fitting range Q 2 < Q 2 lax or Q^ n < Q 2 . This is shown 
in Fig. [6j where we successively limit the fit to smaller (larger) momenta. Note that the 
increasing errors to the left or the right are due to the decrease in the number of fitted data 
points. Within these errors, the change in the monopole mass is consistent with statistical 
fluctuations. From Figs. [5j and we can conclude that the monopole ansatz works well in 
the entire region for which we have lattice data, from Q 2 = to about 4 GeV 2 . 

The results discussed so far have used the lattice data normalised as in ([T9|) . Using 

Z y Fi at < barc (0) = i^ at < ron (0) = 1, (24) 

we can determine Zy from our (unrenormaliscd) data at zero momentum transfer. We find 
reasonable agreement with the values of Zy given in |llj , albeit with errors that are larger 
by at least an order of magnitude. The bigger errors are likely due to our choice of i s ink, 
which results in noisier two-point functions. 
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Figure 5: Effective monopole masses M c ff(Q 2 ) defined in (|23]l . together with the correspond- 
ing monopole masses from Tabie [3] (dotted lines) for a sample of our lattices from small to 
large pion masses (lattices number 8, 7, 4, 9, 2, 10, 15, and 11 from top to bottom). For 
better visibility we omitted 2 points with very large errors in the plot, but included them in 
the fit. 
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Figure 6: Combined fits to (f2"Tj) with reduced fitting ranges in Q 2 . For the left plot Q^ ax is 
decreased, while Qn lin is increased for the right plot. We use bins of 50MeV 2 and show only 
points where the number of data points in the fit of F„ changed. 
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Table 5: Overview of our finite size runs. Note that we use the pion mass and lattice spacing 
of the largest lattice also for the smaller ones. They are given in Table [T] and not repeated 
here. 






K 


# 


N 3 x T 


L [fax] 




M lat [GeV] 


AAf lat [McV] 


5.29 


0.13550 


10 


24 3 x 48 


2.0 


7.8 


0.971(14) 


1.4 






10a 


16 3 x 32 


1.3 


5.2 


0.928(16) 


19.7 






10b 


12 3 x 32 


1.0 


3.9 


0.841(48) 


75.0 


5.29 


0.13590 


11 


24 3 x 48 


1.9 


5.7 


0.854(15) 


12.6 






11a 


16 3 x 32 


1.3 


3.8 


0.786(18) 


90.3 






lib 


12 3 x 32 


1.0 


2.9 


0.513(31) 


263.1 



5.2 Finite volume and discretisation effects 

Let us now turn to the discussion of lattice artifacts. Apart from the extrapolation to the 
physical pion mass there are two more limits to be taken: the infinite volume limit and the 
continuum limit. The large number of lattices available allows us to investigate both. In order 
to study the volume dependence of our results, we make use of two sets of configurations that 
have the same parameters j3, k for the lattice action but different volumes (see Table ■ In 
Fig. [7^, we show the monopolc masses fitted according to Eq. (|2T)|) as a function of the lattice 
size L. We use the pion mass ra, and lattice spacing a determined for the lattice with the 
largest volume also for the smaller ones. Figure [TJd gives an overview of our lattices in the 
m n -L plane. 

To obtain some understanding of the volume dependence one may have recourse to chiral 
perturbation theory. The volume dependence of the pion charge radius has been investigated 
to one-loop order in various approaches of chiral perturbation theory [23 [5^1 • hi the 
continuum limit, the result of the lattice regularised calculation in [27j amounts to a finite 
size correction of 

(^L-^^iE^hD, (25) 

where the sum runs over all three- vectors n ^ with integer components and f n w 92 MeV 
is the pion decay constant. Note that the finite size correction of the charge radius is not 
proportional to m 2 , unlike for other quantities such as the pion decay constant or the nucleon 
axial coupling. The leading contribution in Eq. (|25p for large values of m^L is proportional 
to Kofa-x-L) ~ y7r/ (2m w L) e~ m ^ L , Unfortunately we cannot expect chiral perturbation 
theory to be applicable at the pion masses and lattice volumes used in our simulations. This 
includes the result ([23]) . which we take however as a guide for the functional form of the 
volume dependence. We thus change the monopole mass in ([2Tj) to^l 

M 2 (ml,L) = c Q + Cl mi + C2C- m - L . (26) 

We then perform a combined fit to the data of all lattices in Table [T] except for number 12 
(see below), including in addition the 16 3 x 32 lattices of the finite volume runs (numbers 
10a and 11a). The result is represented by the solid lines in Fig. [7^. The fitted parameters 
are c = 0.553(29) GeV 2 , c x = 0.612(35) and c 2 = -6.97(1.71) GeV 2 at x 2 /d.o.f. = 0.62 

''Taking the Bessel function Ko(m w L) instead of e~~ m7rL does not change our results significantly. 
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Figure 7: (a) Monopole mass vs. lattice size in our finite volume data sets with (3 = 5.29 and 
K = 0.1355 (upper points) or k = 0.1359 (lower points). The curves correspond to a fit to 
(|2"6"|) as discussed in the text, (b) Overview of pion masses and lattice sizes for our complete 
data set. The dotted lines mark our finite size runs. 



which gives M p h ys = 0.751(19) GeV for the infinite volume limit of the monopole mass at 
the physical point. Compared with the value 0.727(16) GeV obtained in the fit (f2"Tj) without 
volume dependence, this represents a small overall finite-size effect. The fitted parameters do 
not change significantly if we only fit the 16 3 x 32 and 24 3 x 48 data sets of the finite volume 
runs, i.e. the data corresponding to the four rightmost points in Fig. [7K (lattices number 10, 
10a, 11, and 11a). We have not included the 12 3 x 32 lattices in the fit (f2"6"| since we cannot 
expect our simple ansatz to hold down to lattice sizes of 1 fm. Qualitatively, our fit is not 
too bad even in this region, as shown by the dotted lines in Fig. [7^. 

With the fitted parameters we can estimate the finite volume shift for each of our lattices 
as given in Table [3] Except for a few lattices we find very small effects. We do not expect 
that with the simple form (|2"6"|) fitted to our finite volume data at — 591 MeV and 
m.^ = 769 MeV (the dotted lines in Fig. [7)d) we can estimate volume effects for pion masses 
as low as 400 MeV. We therefore have excluded lattice number 12 from our finite volume 
investigation. 

Before discussing the scaling behaviour, let us briefly discuss the possibility of 0(a) 
improving the local vector current. The improved current has the form 

V™v{x) = u{x)^u(x) + c v ad v T„ v {x) , 
T^ix) = iu(x)a^ v u{x) . 

The improvement coefficient cy is only known from lattice perturbation theory [10] because 
the only non-perturbative calculations to date are for quenched fermions (see e.g. [2"5]). 
However, even with tadpole improvement the perturbative value for our coarsest lattice is 
cy « —0.027. This is so small that we expect no sizable effect on our results. To see this, we 
plot in Fig. [8] the ratio 

"imp W ) — i , n ,_ / u 
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Figure 8: The ratio ?*i mp (Q 2 ) defined in Eq. (f28|) . evaluated for our coarsest lattice (ft = 5.20, 
K = 0.1342). To obtain the effect of 0{a) improving the current, this ratio needs to be 
multiplied with cy. 

of the pion matrix elements for the two operators on the r.h.s. of Eq. (|27|) . The dependence 
on the index [i cancels in this ratio. Note that here we use unrenormalised lattice data and 
that we still have to multiply with cy in order to obtain the effect of the improvement term 
in the current. This example plot is for our coarsest lattice (/3 = 5.20 and k = 0.1342), where 
the improvement term should have the largest impact. To gain a feeling for the possible size 
of the effect, we used a fixed value of cy = —0.3 to compute the effect on a sub-set of our 
lattices (lattices number 2, 6, 11, 15). Although this improvement coefficient is more than 
ten times larger than the tadpole improved value for our coarsest lattice, the shift of the 
monopole mass was moderate with 6 to 10%. Given the size of our statistical errors on F„ 
and the fact that a reliable value for cy is not known for our lattices, we decided to neglect 
operator improvement and use the local vector current. 

We now investigate the scaling behaviour by extrapolating our values for the monopole 
mass to the physical pion mass separately for each {3 (see the upper plots in Fig.[!|]). We again 
assume a linear relation between the squared monopole and pion masses. The extrapolated 
values can then be studied as a function of the lattice spacing a, using ro/a extrapolated to 
the chiral but not to the continuum limit [H]EI This is shown in the lower plot in Fig. [U 
While the three rightmost data points in the lower plot of Fig. [9] strongly suggest that no 
discretisation errors are present within statistical errors, it requires additional simulation 
points to see if the leftmost data point in the lower plot of Fig. [5] represents a downwards 
trend or is just an outlier. From the discussion above and the overview in Table [3] we recall 
that some of the points at low pion mass may be affected by finite volume corrections. We 
have repeated the fits shown in Fig. with squared monopole masses shifted upwards by 
c 2 e~ m7 ' L , where C2 (in units of r^ 2 ) was taken from the global fit described after Eq. (|26|) . 
Note that the pion mass of 400 McV is excluded from this global fit for the reasons given 
above. The result shows an increase of M p h y s mainly for f3 = 5.20 and 5.40 but is again 

4 We have updated values for ro/a w.r.t. 1291 : for f3 = 5.20, 5.25, 5.29, and 5.40 they are ro/a = 
5.444(72), 5.851(85), 6.158(53), and 6.951(54), respectively. 
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Figure 9: Scaling test: the upper plots show extrapolations as in Fig.[3]for each (3 separately. 
The lower plot shows the extrapolated values of the monopole mass at the corresponding 
square of the lattice spacing. 



consistent with no a dependence. Given the lever arm in o 2 and the size of our statistical 
and finite size errors, we refrain from including an explicit a dependence of the monopole 
mass in our global fit (l2"Tj) . 



6 Conclusion 

We have calculated the electromagnetic form factor of the pion, using lattice configurations 
generated by the QCDSF/UKQCD collaboration with two flavours of dynamical, 0(a) im- 
proved Wilson fermions. The corresponding pion masses range from 400 to 1180 MeV. The 
momentum dependence of the pion form factor was studied up to Q 2 around 4 GcV 2 . Within 
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Table 6: An overview of lattice results for the pion charge radius along with the experimental 
value. We only quote results that are extrapolated to the physical point. The quoted lattice 
errors are purely statistical. 



(r 2 ) [fm 2 ] 


type of result 


Reference 


0.452(11) 


experimental value 


PDG 2004 [21] 


0.441(19) 
0.396(10) 
0.37(2) 
0.310(46) 


Clover improved Wilson fermions, Nt = 2 
Clover improved Wilson fermions, Nf = 2 
Wilson fermions. quenched 
hybrid ASQTAD/DWF, N f = 2 + 1,3 


this work 
JLQCD [5] 



LHPC [E] 



errors, the pion form factor is described very well by a monopole form ([2H|) in this range, 
for all our lattice pion masses. A linear chiral extrapolation to the physical pion mass leads 
to a monopole mass of M = 0.727(16) GeV. This corresponds to a squared charge radius 
(r 2 ) = 0.441(19) fm 2 , in good agreement with experiment. Our extrapolated lattice data for 
the form factor is compared with experimental measurements in Fig. fj] Other lattice results 
are quoted in Table 

The large parameter space of the gauge configurations we used makes it possible to explore 
artifacts arising from the finite lattice spacing and volume. An empirical fit allowing for a 
volume dependence leads to an increase of the monopole mass by 3% at infinite volume and 
the physical point. Within errors, our results show no clear dependence on the lattice spacing 
in the range a = 0.07 - 0.11 fm of our simulations. Including estimates for systematical 
errors, our final result then is M = 0.727 ±0.016 (stat)± 0.046 (syst) + 0.024 (vol) GeV, which 
translates to a charge radius of (r 2 ) = 0.441 ± 0.019 (stat) ± 0.056 (syst) - 0.029 (vol) fm 2 . 
The first error is purely statistical, followed by a systematic uncertainty due to the ansatz 
for the fitting function and the extrapolation to physical pion masses (for which we added 
in quadrature the errors AM GXt and AMfi t obtained in Section [5~Tj) . The last error reflects a 
possible shift because of finite volume effects as just discussed. We have set the scale using 
the Sommer parameter with ro = 0.467 fm. We note that the analysis leading to our result 
for M is independent of the scale setting, so that a different value of ro would lead to a simple 
rescaling of the above values. 
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